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Abstract 

We address a number of aspects of the arrival time problem defined using a complex potential of 
step function form. We concentrate on the limit of a weak potential, in which the resulting arrival 
i-£h ' time distribution function is closely related to the quantum- mechanical current. We first consider 



> 

On 



ON 

o 



the analagous classical arrival time problem involving an absorbing potential, and this sheds some 
light on certain aspects of the quantum case. In the quantum case, we review the path decomposi- 
tion expansion (PDX), in which the propagator is factored across a surface of constant time, so is 
very useful for potentials of step function form. We use the PDX to derive the usual scattering wave 



functions and the arrival time distribution function. This method gives a direct and geometrically 
appealing account of known results (but also points the way to how they can be extended to more 



g ; genera! complex potentials). We n s e these results to carry out a decoherent hrstor.es analys.s of 



the arrival time problem, taking advantage of a recently demonstrated connection between pulsed 
measurements and complex potentials. We obtain very simple and plausible expressions for the 
class operators (describing the amplitudes for crossing the origin during intervals of time) and 
show that decoherence of histories is obtained for a wide class of initial states (such as simple wave 
packets and superpositions of wave packets). We find that the decoherent histories approach gives 
results with a sensible classical limit that are fully compatible with standard results on the arrival 
time problem. We also find some interesting connections between backflow and decoherence. 
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I. INTRODUCTION 



A. The Arrival Time Problem 

The arrival time problem has attracted some considerable interest in recent years In 
the one-dimensional statement of this problem, one considers an initial wave function \ip) 
concentrated in the region x > and consisting entirely of negative momenta. The question 
is then to find the probability U(t)cIt that the particle crosses x = between time r and 
r + dr. (See Figure 1.) 

Two particular candidate expressions for the arrival time distribution II (r) are central to 
the discussion. First, there is the current density 

where IV'/C 7 ")) ^ s the freely evolved state, which arises from elementary considerations of 
the Schrodinger equation. (For convenience we work in units in which H = 1). This is 
sensible classically but can be negative in the quantum case for certain states consisting of 
superpositions of different momenta . Second, a simple operator re-ordering of J(t) 

gives the "ideal" arrival time distribution of Kijowski 



n*(r; 

m 



i^/(r)||p| 1/a ^)|p| 1/a |^(r)> (1.2) 



which is clearly positive. Both of these distributions are measurable 



a. 



B. Complex Potentials 

An interesting question is the extent to which such expressions emerge from more elab- 
orate measurement or axiomatic schemes. There are many such schemes. Here, we will 
focus on expressions for the arrival time distribution arising from the inclusion of a complex 
potential 

V(x) = -iV Q 9(-x) (1.3) 
in the Schrodinger equation. With such a potential, the state at time r is 

|V(r)) = exp (-iH T - V 6(-x)t) |^> (1.4) 
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where Hq is the free Hamiltonian. The idea here is that the part of the wave packet that 
reaches the origin during the time interval [0, r] should be absorbed, so that 

JV(r) = (V(rMr)) (1-5) 

is the probability of not crossing x = during the time interval [0, r]. The probability of 
crossing between r and r + dr is then 

Complex potentials such as Eq. fll.3l) were originally considered by Allcock in his seminal 
work on arrival time [8J and have subsequently appeared in detector models 6, 9]. (See also 
Ref. for further work with complex potentials). A recent interesting result of Echanobe 

et al. is that under certain conditions, evolution according to Eq. (jl.4p is essentially the same 
as pulsed measurements, in which the wave function is projected onto x > at discrete time 



intervals 



For large Vo, the wave function defined by the evolution Eq. fll.4p undergoes significant 
reflection, with total reflection in the "Zeno limit" , Vq — > oo [l2j . Here, we are interested in 
the opposite case of small Vo, where there is small reflection and Eq. fll.6p can give reasonable 
expressions for the arrival time distribution. A number of different authors p, Q, [§| indicate 
that the resulting distribution is of the form 



n(r) = / dt R(V ,T-t)J(t) (1.7) 

J — oo 

where J(t) is the current Eq. fll.ip and 

R(V Q , t) = 2V 6(t) exp (-2V t) (1.8) 

It is therefore closely related to the current J(t) but also includes the influence of the complex 
potential via the "apparatus resolution function" R(V ,t). 

The first aim of this paper is to look in some detail at the calculation and properties of the 
arrival time distribution Eq. (11.71) defined using a complex potential. In particular, we will 
use path integral methods, which in some ways give are more concise and transparent than 
previous derivations (and also suggest generalizations to complex potentials more general 
than Eq. fll.3p ). We will also explore some of the properties of the result, Eq. fll.7p . 

The second aim of this paper is to carry out a decoherent histories analysis of the arrival 
time problem. This turns out to be closely related to the complex potentials definition of 



arrival time and was in fact the original motivation for exploring complex potentials. In 
brief, our aim is to see if standard results for n(r), such as Eq. fll.ip or Eq. (11.71) have a 
natural place in the decoherent histories approach. 



C. The Decoherent Histories Approach to the Arrival Time Problem 



We now briefly review the decoherent histories approach to the arrival time pro blem 



Q Q, Q. 
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18} . In the decoherent histories approach 
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probabilities are assigned to histories via the formula, 

p(a u a 2 , ■ ■ ■ ) = Tr (C„pC^) 



;i.9) 



where the class operator Ca denotes a time-ordered string of projectors P a interspersed with 
unitary evolution, 



r< p „— iH(tn— tn-l) p a—iH(t2—ti) p 



1.10) 



and a denotes the string a±, a 2 , ■ • • a n . The class operator satisfies the condition 

J2C* = e~ iHT (1.11) 



where r is the total time interval, r = t n — t\. Interference between pairs of histories is 
measured by the decoherence functional, 



It satisfies the relations, 



D(a,d) = D*(a',a) 
yj D{a,a) = 1 



a,a' 



^D(a,a) = ^2p(a) = 1 



(1.12) 

(1.13) 
(1.14) 

(1.15) 



Of particular interest are sets of histories which satisfy the condition of decoherence, which 
is 

%,a') = if a^a' (1.16) 

Decoherence implies the weaker condition of consistency, which is that KeD(a, a') = for 
a 7^ a', and this is equivalent to the requirement that the above probabilities satisfy the 



probability sum rules. In most situations of physical interest, both the real and imaginary 
parts of D(a,c/) vanish for a ^ a'ja condition we shall simply call decoherence, and is 



related to the existence of records 
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251 ]. Decoherence is only approximate in general 



which raises the question of how to measure approximate decoherence. The decoherence 



functional satisfies the inequality 24] 



]£>(", a') I 2 <p{oL)P{d) (1.17) 
This suggests that a sensible measure of approximate decoherence is 

\D(a,g/)\ 2 < p(o)j)(a') (1.18) 
The approach also permits other types of class operators which are not simply strings of 



projections, but sums of such strings. These are often called inhomogeneous histories 26] 
and are relevant to questions involving time in a non-trivial way. For example, for a given 
class operator the object 1 — Ca is also a class operator but it not equal to a simple 
string of projections. Unlike homogenous histories, inhomogeneous histories do not satisfy 
condition Eq. fll.15l) in general, except when there is decoherence. 

A quantity closely related to the probabilities is the quasi-probability 

q(a) = Tr (C«pe iHT ) (1.19) 

Using Eq. fll.lll) . this satisfies 

= p(ql) + D (^^) (i- 20 ) 

(where a is fixed in the sum). This means that when there is decoherence, the probabilities 
for histories p(a) are equal to the simpler expression q(a). (The converse is generally not 
true except for in very simple circumstances). Note that q(a) is not positive (or even real) 
in general, but it is positive and real when there is decoherence. These facts turn out to be 
relevant to our analysis of the arrival time problem. 

We ultimately seek a decoherent histories account of the arrival time probability U(t)dt, 
the probability for the particle to arrive in an infinitesimal interval [t,t + dt]. However, for 
simplicity, we first consider the simpler problem of computing the probability for arriving 



Some earlier papers on the deco 
of the class operators 14 



during a finite (possibly large) interval [0, r]. We consider an incoming state entirely localized 
in x > 0, and we partition the system's histories into two classes: histories that either cross 
or do not cross x = during the time interval [0, r]. We seek class operators C c , C nc 
corresponding to these two classes of histories. 

herent histories approach adopted the following definition 
15! . If] . \\\ . (This definition is problematic, as we shall see, 
but sets the stage for the corrected version we shall use here.) Let P = 9(x) denote the 
projection onto the positive x-axis. To define the class operator for histories which do not 
cross x = 0, we split the time interval into N parts of size e, and the class operator is defined 
by 

C nc = \imPe- iHe P---e- iHt P (1.21) 

where there are N + 1 projections and N unitary evolution operators, and iV — > 00 in such 
a way that r = iVe is constant. The limit actually yields the so-called restricted propagator, 

C nc = g r (r,0) (1.22) 

This object is also given by the path integral expression 

{xi\g r (r, Q)\x ) = / Vxexp(iS) (1-23) 



where the integral is over all paths from x(0) = Xq to x{r) = x\ that always remain in 
x(t) > 0. The class operator for crossing the surface is then defined by 



C c = e~ lHT - C nc (1.24) 

However, as indicated, there is a problem with this definition. The class operator 
Eq. ( 11.211 ) suffers from the quantum Zeno effect 12j - projecting continually in time onto the 
region x > prevents the system from leaving it and the probability for not crossing x = 
is unity, 

Pnc = Tr (C ncP Cl) = 1 (1.25) 

for any initial state. This is easily seen from the observation that the restricted propagator 
defined by the limit of Eq. (11.211) may actually be written, 

g r (r, 0) = P exp (-iPHPr) (1.26) 
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so is unitary in the Hilbert space of states with support only in x > [27|, |28|]. Differently 
put, an incoming wave packet evolving according to the restricted propagator undergoes 
total reflection, so never crosses x — 0. These results are clearly unphysical and have no 
sensible classical limit. 

The problem here is that the system is monitored too closely to allow the wave packet 
to actually cross x = 0. The resolution is therefore to relax the monitoring in such a way 
that something interesting can happen. There are two obvious ways in which this may be 
achieved. 

The first option is to simply decline to take the limit e — > in Eq.f ll.21i ). so we define the 
class operator for not crossing to be, 

C e nc = Pe~ iHe P ■ ■ ■ e~ iHe P (1.27) 

where again r = iVe, but N is finite and e > 0. Clearly if e is large enough the system will 
be monitored sufficiently infrequently to let the wave packet cross x = without too much 
reflection. Studies of the quantum Zeno effect suggest that the appropriate lower limit on e 
is 



where H n is the free Hamiltonian 



£ > (L28) 



The second option is to retain the limit in Eq. fll.21j) . but "soften" the projections to 
POVMs, that is, to replace P = 8(x) with a function which is approximately 1 for large 
positive x, approximately for large negative x, and with a smooth transition in between. 

What we will do in this paper is adopt a definition of the class operator which involves ele- 
ments of both of these options. In particular, we exploit the useful observation of Echanobe et 



al 



ll| . which is that the string of projections interspersed with unitary evolution, Eq. fll.27p 



is approximately equivalent to evolution in the presence of the complex potential Eq. fll.3j) . 

P e - iHe P---e- iHe P « exp (-iH Q T - V q 9(-x)t) (1.29) 
This connection follows from first noting that 

p = e(x) w exp (-V o 0(-x)e) (1.30) 

as long as 

V e > 1 (1.31) 



We now note that the approximate equivalence 

exp {-iH e) exp (-V 9(-x)e) « exp {-iH e - V 9(-x)e) (1.32) 
will hold as long as 

V e 2 \([H ,9(-x)])\^^ (1.33) 
Echanobe et al. put an upper bound on the left-hand side using the Schrodinger-Robertson 
inequality and the two inequalities together read 

1 « V 2 e 2 « (1.34) 

which can be satisfied as long as 

V > AH (1.35) 

(Note that the Zeno limit restriction Eq. fll.28p is not in fact relevant here since the quantum 
Zeno effect goes away when exact projectors are replaced by quasi-projectors, as in Eq. fll.30p 

H). 

With a pinch of salt, these arguments indicate that a connection something like Eq. fll.3ip 
will hold. This equivalence and demonstration of it is somewhat heuristic, but appears to 



be backed up by numerical evidence 
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29] . We will accept it in what follows, but it would 
be of interest to obtain a more rigorous proof of this connection. 

This is a very useful result, since the problem of evolution in the presence of the complex 
potential Eq. (11.31) is straightforward to solve, but the evolution described by Eq. fll.27p could 
be difficult to solve analytically. In our decoherent histories analysis of the system, we will 
therefore define the class operator for not crossing to be 

C nc = exp {-iH r - V 9{-x)r) (1.36) 

with the crossing class operator defined by Eq. (11.241) . These definitions will be extended to 
class operators for crossing during any one of a set of small intervals of size A. 

These are the class operators we need to give a decoherent histories account of the origin of 

n(r). 

Our main result is that for intervals of size A ^> 1/Vo, these class operators are given 
approximately by 

ptk+i (—1) 



-lH T 



(9(x(t k )) - 9(x(t k+ i))) (1.37) 



Significantly, the dependence on the complex potential has dropped out entirely. We will 
show that there is decoherence for an interesting class of states, and, for such states the 
probabilities are then given by Eq. (11.191) . which has the form 



fe+i 



q(t k ,t k+ i) = / dt J(t) (1.38) 



This, we will show, coincides with the expected result Eq. (11.71) . when is integrated 
over a range of time much greater than 1/Vq. Hence there is complete agreement with 
standard results on the arrival time distribution at sufficiently coarse-grained timescales. 
Furthermore, our results shed some light on the problem of backflow - we find that the 
situations when Eq. (11.381) is negative are those in which there is no decoherence in which 
case probabilities cannot be assigned. 

In another paper, we compute the crossing class operators in a simpler but more heuristic 
way, by exploring a semiclassical approximation to Eq. fll.27p keeping e finite. 30j. The 
results are essentially the same. 



D. This Paper 

The initial motive and overall aim of this paper are to give a decoherent histories account 
of the arrival time question which has a sensible classical limit. As stated, we seek a deco- 
herent histories analysis account of standard results, such as Eq.f ll.7p . However, along the 
way we have found some new derivations of known results and will explore these areas too. 

In Section 2, we consider the classical analogue of the arrival time defined by a complex 
potential. This sheds light on the form of the result Eq. f ll .71) and in particular, the origin of 
the resolution function Eq.f ll.8p . In Section 3, we review the path decomposition expansion 
(PDX), which is a useful path integral technique for factoring the propagator across the 
surface x = 0, so is very useful for systems with a step function potential. In Section 4, 
we use the PDX to derive the scattering wave functions corresponding to evolution with a 
complex step potential. These are of course known, but the PDX gives a useful and concise 
derivation of them. In Section 5, we use the PDX to calculation the arrival time distribution 
function Eq. fll.6p . and very readily obtain the expected result of the form Eq.f ll.7p . In 
Section 6 and 7, we use the results of the previous sections to carry out the decoherent 
histories analysis, as outlined in the previous subsection. We summarize and conclude in 
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Section 8. 



II. THE CLASSICAL CASE 

Before looking at the quantum arrival time problem, it is enlightening to look at the 
corresponding classical arrival time problem defined using an absorbing potential. This 
gives some understanding of the expected form of the result in the quantum case, Eq. (11.71) . 

We consider a classical phase space distribution w t (p,q), with initial value wo(p,q) con- 
centrated entirely in q > with only negative momenta. The appropriate evolution equation 
corresponding to the quantum case Eq. (11.41) is 

— = -JL 2 V(q)w (2.1) 

at m oq 

where V(q) = Vq8(— q). This form may be deduced, for example, by computing the evolution 
equation of the Wigner function corresponding to Eq. (ll.4p and dropping the higher order 
terms (involving powers of K). Eq. (12 . 1 p is readily solved and has solution 



w t{Pi q) — exp ^— 2 J ds V(q — ps/m) J wo{p,q — pr /m) (2.2) 
Following the corresponding steps in the quantum case, the survival probability is 

/oo />oo 
dp dqw T (p,q) (2.3) 
-oo J —oo 



and the arrival time distribution is 



dN 



o 



n(r) = - — = 2V Q / dp dq w T (p, q) (2.4) 

J —oo J —oo 

where we have made use of Eq. ( 12. ID . This expression is conveniently rewritten by noting 
that, again using Eq. (I2.1[) . Il(r) obeys the equation 

dll f°° v 

— + 2V U = -2V / dp ^w T {p, 0) (2.5) 

(XT J -oo m 

This may be solved to yield 

n(r) = -2V Q f dt e" 2 ^^"*) f dp —w t (p, 0) (2.6) 
Jo J-oo rn 

From Eq. (12.21) . we see that 



w t{Pi 0) = exp ( — 2Vo / ds 6{ps/m) J wo(p, —pt/m) (2.7) 
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but since the momenta are all negative the exponential factor makes no contribution. We 
thus obtain 



n(r) = 2V [ T dt e - 2V °^-') J(t) (2.8) 
Jo 



where 



J(t) = - [ dp ^w (p, -pt/m) (2.9) 

The current J(t) is the usual classical arrival time distribution that would have been obtained 
in the absence of the absorbing potential. 

The result Eq. fl2.8D is very close in form to the expected quantum result, Eq. fll.7l) . differing 
only in the range of integration. The lower limit of — oo in Eq. (11.71) arises as a result of the 
approximations used in deriving it and this difference is not significant since we expect the 
current to be approximately zero anyway for t < 0. 

The classical result Eq.f l2.8p helps to understand the role of the resolution function R(t) 
in the quantum case Eq. fll.7p - it in essence describes the classical influence of the absorbing 
potential used to model the detector. More precisely, R(t) is actually related to a sort of 
coarse graining in time, and this we now demonstrate. Eq. fl2.8l) gives the probability n(r)<ir 
for arriving during the infinitesimal time interval [r, r + dr]. Supposer we consider the 
probability for arriving during a finite time interval, [n, r 2 ]. This is given by 

p(t 2 ,ti) = dr n(r) 

Jti 

= f 2 dr [ dt 2V e - 2V!, (--*) J(t) (2.10) 
Jt! Jo 

Rearranging the order of integration and integrating over r yields, 
p(r 2 , ri) = f'dtfdr 2V e - 2V °^-*) J(t) 

Jo Jti 

+ dt j T ' dr 2V e- 2Vo( - T ~V J(t) (2.11) 

= P dt (e~ 2y "^-') - e -2Vo(^-*)) 
Jo 

+ f 2 dt (1 - e - 2V ^ T2 - t] ) J(t) (2.12) 

We will see in what follows that 1/Vq plays a role as a fundamental short timescale in the 
problem. So now suppose we assume that Ti, r 2 and (r 2 — r{) are all much greater than 
1/Vq. It follows that all the exponential terms may be dropped in Eq. fl2.12l) and we obtain 



11 



the very simple result, 

p(r 2 ,ri) ^ I" dt J{t) (2.13) 

Jti 

That is, all dependence on the resolution function R(t) and the complex potential parameter 
Vo completely drops out when we look at probabilities defined on timescales much greater 
than 1/Vq. This result is very relevant to the decoherent histories analysis considered later 
where it is natural to look at the arrival time during a finite time interval. 



III. THE PATH DECOMPOSITION EXPANSION 



In this section we describe some useful path integral techniques. We wish to evaluate the 
propagator 

g(x ± , t\x , 0) = exp (-iH r - V d(-x)f(x)r) \x ) (3.1) 
for arbitrary x\ and Xq > 0. This may be calculated using a sum over paths, 

g(xi, t\xq, 0) = J Vxexp(iS) (3.2) 

where 

S = jf dt Qmx 2 + iV Q e{-x)f\x)^j (3.3) 

and the sum is over all paths x(t) from x(0) = xq to x(t) = x\. 

To deal with the step function form of the potential we need to split off the sections of the 
paths lying entirely in x > or x < 0. The way to do this is to use the path decomposition 



expansion or PDX 3jJ, |32j, [33]. Consider first paths from xq > to x\ < 0. Each path from 



Xq > to X\ < will typically cross x = many times, but all paths have a first crossing, 
at time t\, say. As a consequence of this, it is possible to derive the formula, 

i [ T dg r 
flf(a?i,r|a5 ,0) = — J dt x g(x 1: r|0, ti) — (x, tx\x , 0) \ x=0 (3.4) 

Here, g r (x,t\xo,Q) is the restricted propagator given by a sum over paths of the form (13. 2\i 
but with all paths restricted to x(t) > 0. It vanishes when either end point is the origin but 
its derivative at x = is non-zero (and in fact the derivative of g r corresponds to a sum over 
all paths in x > which end on x = |32j). 

It is also useful to record a PDX formula involving the last crossing time t 2 , for xo > 
and x\ < 0, 

i [ T dg r 

g(x 1 ,r\x o ,0) = -— J dt 2 -q^{ x i, t \ x M)\ x=q #(0, t 2 ko, 0) (3.5) 
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These two formulae may be combined to give a first and last crossing version of the PDX, 

i r dt f'\> 9 3 , , 

'o Jo 



g(xi,r\x ,0) = / dt 2 / dt% —(xi,r\x,t 2 )\ x=0 #(0,t 2 |0,ti) -^-(ar, ti|aco, 0) | x=0 



(3.6) 

This is clearly very useful for a step potential since the propagator is decomposed in terms 
of propagation in x < and in x > 0, essentially reducing the problem to that of computing 
the propagator along x = 0, g(0,t2\0,ti). (See Figure 2.) 

For paths with x > and x% > 0, Eq. fl3.4p is modified by the addition of a term 
g r (xi, t\xo, 0), corresponding to a sum over paths which never cross x = 0, so we have 

g(x 1 ,r\x o ,0) = — y dti #(xi, r|0, ti)-^(x, ti|2 , 0))^ + g r (xi, t|x , 0) (3.7) 

Again a further decomposition involving the last crossing, as in Eq. (13.61) can also be included. 

The various elements of these expressions are easily calculated for a potential of simple 
step function form V(x) = Vq6(— x). The restricted propagator in x > is given by the 
method of images expression 

g r {xx, r\x , 0) = 9{xx)9{xq) (gf(xi, r\x , 0) - gf(-x h r\x , 0)) (3.8) 

where gf denotes the free particle propagator 

g f {x u t\x , 0) = (— j exp ^ j (3.9) 

It follows that 

^{xM\^)\ x=Q = 2^(0,^1x0,0)^X0) (3.10) 

The restricted propagator in x < is given by Eq. fl3.8p . multiplied by exp(— V t). The only 
complicated propagator to calculate is the propagation from the origin to itself along the 



edge of the potential, and in the case V(x) = Vq9(— x) this is given by [341 ]. 



/ , x /m\V2 (1 - exp(-V t)) 

Using these results we may write down the full solution to the evolution with a complex 
potential, described by Eq. (13. II) . for an initial state ip(x) with support only in x > and 
with negative momenta. It has the form, 

i/;(x, r) = 9(-x)il) tr (x, r) + 9{x) (i[>ref(x, r) + i[) f (x, r)) (3.12) 
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Here, iptr is the transmitted wave function and is given by the propagation of the initial 
state ijj{x) using the PDX formulae Eq. fl3.4p or Eq. fl3.6p . 

The remaining part ipref + ipf is the wave function obtained by propagating using the 
PDX formula for initial and final points both in x > 0, Eq. fl3.7l) (rewritten using Eq. fl3.6p if 
necessary). It is appropriate to break this into the two pieces ip re f, ipf defined as follows: the 
reflected wave function ip re f consists of that obtained using the first term in Eq. (I3.7j) together 
with the reflected part — <?/(— X\, t\xq, 0) of the restricted propagator, g r . This definition 
ensures that ip re f — > as the complex potential in x < goes to zero. The remaining part, 
if}f, is the other part of the restricted propagator so is simply free propagation in i > 0. 
This corresponds to the part of the incoming wave packet that never reaches x = during 
the time interval [0, r]. This part clearly goes to zero for large r. 

IV. SOLUTION THROUGH STATIONARY SCATTERING STATES 

In this section we use the PDX to derive the standard representations of the scattering 
solutions to the Schrodinger equation with the simple complex potential Eq. fll.3p . These are 
known results but this derivation confirms the validity of the PDX method and also allows 
a certain heuristic path integral approximation to be tested. The results will also be useful 
for the decoherent histories analysis in Section 6. 

The transmitted and reflected wave functions are defined above in Eq.f l3.12l) . For large 
r, a freely evolving packet moves entirely into x < so that the free part, ipf(x,r) is 
zero, leaving just the transmitted and reflected parts. Following the above definition, the 
transmitted wave function is given by 

1 f T PT—S 

iptr(x,r) = — 2 ds dv (x\ exp (— iH s) p\0) e~ Vos 



where |0) denotes the position eigenstate \x) at x = 0. Also, we have introduced s = t — t± 
and v = t 2 — t±, and H = H — iVo9(—x) is the total Hamiltonian. The scattering wave 
functions concern the regime of large r, so we let the upper limit of the integration ranges 
extend to oo. 

Writing the initial state as a sum of momentum states \p), and introducing E = p 2 /2m, 




x (0| exp (-iHv) |0) (0|pexp {-iH (t - v - a)) 



(4.1) 
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we have 



77T 



ip tr (x,T) = — - dp ds{x\ exp (— iH s) p\0) e 



i(E+iV )s 







x / dv (0\exp(-iHv)\0) e lhv p(0\p) e~ lhT ij(p) 
'o 



To evaluate the s integral, we use the formula [351 ]. 



ds 



m \ V 2 



2-7T2S 



exp % 



As + 



rax 
~2^~ 



— ^ 2 exp (i|x|\/2mA 



(4.2) 



(4.3) 



from which it follows by differentiation with respect to x and setting A = E + iVo that 

ds(x\exp(-iH s)p\0) e i{E+iVo)s = mexp (i\x\[2m(E + iV ))^) (4.4) 



The f integral may be evaluated using the explicit expression for the propagator along the 
edge of the potential, Eq. fl3.lip . together with the formula, 



2iri 



'I - e- Vov ) 



'2m 



(E + iV )*+E2 



We thus obtain the result, 

ip tr (x,r) = 

where 



dp 
v/27 



exp ^— ix[2m(E + iV )] 2 — iEr^j iptr(p) 



~T7 ipip) 



(4.5) 



(4.6) 



(4.7) 



{l + E-^E + tVo) 

Note that in this final result, it is possible to identify the specific effects of the different 
sections of propagation: the propagation along the edge of the potential corresponds to the 
coefficient in the transmission amplitude Eq. (14.71) (which is equal to 1 when Vq = 0), and 
the propagation from final crossing to the final point produces the Vq dependence of the 
exponent. These observations will be useful below. 

The reflected wave function ip re f is defined above using the PDX Eq. fl3.7p (rewritten using 
Eq. (13.61) . The first term in Eq. fl3.7p . the crossing part, is the same as the transmitted case, 
Eq. (14.21) . except that Vq = in the last segment of propagation, from x = to the final 
point, and also the sign of x is reversed. We must also add the effects of the reflection part 
of the restricted propagator, and this simply subtracts the reflection of the incoming wave 
packet. The reflected wave function is therefore given by 

dp 



Aef(x, r) 



exp(ixp — iEr) 4> re f(p) 



(4- 
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where 



Aefip) = Mp) - ip{p) 

l-E-v(E + iV )*) 
(l + E-h(E + iV )^) 



(4.9) 



We thus see that the PDX very readily gives the standard stationary wave functions |8( , 
without having to use the usual (somewhat cumbersome) technique of matching eigenfunc- 
tions at x = 0. In fact, this procedure is in some sense already encoded in the PDX. 

We now use these exact results to check the validity of a useful approximation in the 
path integral representation of the propagator. In the PDX, Eq. fl3.4l) . the awkward part to 
calculate (especially for more general potentials) is the propagation from x = to the final 
point x\ < 0. The exact propagator for this section consists of propagation along the edge 
of the potential followed by restricted propagation from x — to xi, as used in Eq.f l3.6l) . 
However, for sufficiently small Vo, one might expect that in the path integral representation 
of the propagator, the dominant contribution will come from paths in the neighbourhood of 
the straight line path from x = to x < 0. These paths lie almost entirely in a; < 0, and 
one might expect that the propagator is therefore given approximately by 



It is not entirely clear that this is the case, however. On the one hand, the usual semiclassical 
approximation indicates that paths close to the straight line paths dominate, but on the 
other hand, paths in x < are suppressed, so maybe the wiggly paths that spend less time 
in x < make a significant contribution. Since this approximation is potentially a useful 
one, it is useful to compare with the exact result for the transmitted wave packet calculated 
above. 

We therefore evaluate the following approximate expression for the transmitted wave 
function, 



This is the PDX, Eq. fl3.4l) . in operator form with the semiclassical approximation described 
above and we have set s — r — 1\. We now take r — > oo in the integration and evaluate. 
The key integral is, 




(4.10) 




(4.11) 



/ 

Jo 



■oo 



ds (x\e 



—iHos 




2(E + iV ) 



rn 



exp (-ix[2m(E + iV )}^) (4.12) 
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where we have used Eg. ( 14.31) (and recall that x < 0). The resulting expression for the 
transmitted wave function is of the form Eq. fl4.6p . with 

MP) = ~ i,J .,^M P) (4-13) 

This agrees with the exact expression for the transmission coefficient Eq. fl4.7p only when Vq = 
0, with the difference of order Vq/E for small Vq. This establishes that the approximation 
is valid for Vq much less than the energy scale of the initial state. 

V. CALCULATION OF THE ARRIVAL TIME DISTRIBUTION 

With the general complex potential Eq. (11.31) . the arrival time distribution Eq. fll.6p is 
given by 

n(<r) = 2V o (Vv|0(-z)|Vv> (5.1) 

where 

|Vv) = exp (-Mt) 1*0) 

= exp (-iH T - V 9(-x)t) |V> (5.2) 

(so we use H = H — iVq9(— x) to denote the total non-Hermitian Hamiltonian) . We are 
interested in calculating this expression for the case in which Vq is much smaller than the 



energy scale of the initia 
been explored elsewhere 



. state. (The very different limit, of Vq — > oo, the Zeno limit, has 
36|]). 

One way to evaluate Eq.f l5.ip is to use the transmitted wave functions, Eq. (14.21) . However, 
we give here instead a different and more direct derivation using the PDX. We use the first 
crossing PDX, Eq. (13.41) . which is conveniently rewritten as the operator expression, 

1 r 

{x\ exp(— Mt)\^) = / dt (x\ exp(— iH(r — t)) S(x)p exp (— iH t) (5.3) 

m Jo 

Now note that the operator 5(x) = |0)(0| has the simple property that for any operator A 

5(x)A5(x) =5(x)(0\A\Q) (5.4) 

(where, recall, |0) denotes the position eigenstate \x) at x = 0). Inserting Eq. fl5.3p in 
Eq. fl5.ll) . together with the property Eq. fl5.4p and the change of variables s = r—t, s' = r—t', 
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yields 



m 2 



o 

JO J -oo 



nfr) = ^4 I ds' I ds I dx 



x (0| exp (iifV) \x)(x\exp (—iHs) |0) 

x exp (iH (r — s')) p 5(x) p exp (— iH (r — s)) \ip) (5.5) 

We are aiming to show that this coincides with Eq. fll.7p with the current Eq. fll.ip . and the 
main challenge is to show how the p5(x)p combination turns into the combination p5(x) + 
5(x)p in the current Eq. fll.ip . 

Consider first the x integral. Since we are assuming small Vq, we may use the semiclassical 
approximation Eq. fl4.10p . which reads 

/ m \ l/^ ( hit 2 \ 

(x\ e xp(-tHs)\0) « (— ) expU^-Vos) (5.6) 

The x integral may now be carried out, with the result, 

™ 2 k Jo ^m) ( S - S ')i 
x (ip T \ exp (—iH s')p S(x) p exp (iH s) \ip T ) (5.7) 

where denotes the free evolution of the initial state. 
We now carry out one of the time integrals. Note that, 

PT PT PT PT PT 

ds' \ ds= ds' ds+ ds ds' (5.8) 

Jo Jo Js' Jo Js 

In the first integral, we set u = s', v = s — s', and in the second integral we set u — s, 
v — s' — s. We thus obtain 



x 



— (Vv| exp (-iH u)p S(x) p exp (iH (u + v)) \ip T ) 



+ — —{i'A exp (—iH (u + v))p 8{x) p exp (iH u) \ip T ) 



(5.9) 



H) 5 

The factors of l/(±i)a in front of each term come from a careful consideration of the square 
root in the free propagator prefactor (and must have this form because II (r) is real). 

We will assume that Vqt 3> 1, which means that the integrals are concentrated around 
u — v — 0. This means that we may take the upper limit of the v integral to be oo, and it 
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may be carried out, to yield, 



U(t) = 2V n / du e 



-2V u 







2m 

where the operator £(p) is given by 



x ^-(^Pr-u\ V S(x) E(p) + Et(p) p |W-„) (5.10) 



[2m(i/ + 2Vo)]2 



For Vo much less than the energy scale, 

S(P) - || (5-12) 

so S(p) is simply the sign function of the momentum, which is —1 in this case, since the 
initial state consists entirely of negative momenta. Finally, writing u = t — t, we obtain 

n(r) = 2V [ T dt e- 2Vo(r -^ t^-(,/; t \{pS(x) + 5{x) p) \^ t ) 

2m 



o 



T 



2V Q / dt e- 2 ^-^ J(t) (5.13) 



o 



We therefore have precise confirmation of the classical result Eq.f l2.8p . and also agreement 
with the expected quantum result, Eq. fll.7l) . modulo the issues already discussed concerning 
the range of integration of t. 

Some comments are in order concerning the positivity of the result for II (r). The ex- 
pression (15.101) is positive because it was derived from the manifestly positive expression 
Eq. ( 15.51) . (Two approximations were used: the semiclassical approximation Eq. (14.101) . and 
the condition Vqt ^> 1, neither of which affect the positivity of the result.) 

However, to obtain the final result Eq. ( 15.131) we took the limit Vo — > in the current part 
only of Eq. (15. 101) . leaving behind the Vo-dependent term in the exponential part, and the 
resulting expression is not guaranteed to be positive. In Eq. (l5.13[ ). J(t) is not always positive 
due to the backflow effect 2|, |3] and integration over time does not necessarily remedy the 
situation. (See Appendix A for a more thorough discussion of this point). The lack of 
positivity for a II (r) obtained in this way is not surprising since taking the limit Vo — > in 
one part of the expression Eq. (15.10j) only but not the other will not necessarily preserve its 
property of positivity. The violation of positivity is generally small, however, so Eq. (l5.13p 
may still be a good approximation to the manifestly positive expression Eq. (l5.10p . 
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It should also be added that it would be misleading to explore the first order corrections 
in Vo in going from Eq. fl5.10p to Eq. fl5.13l) . since comparable correction terms have already 
been dropped in using the semiclassical approximation Eq. fl4.10p . 



VI. DECOHERENT HISTORIES ANALYSIS FOR A SINGLE LARGE TIME IN- 
TERVAL 

We now consider the decoherent histories analysis of this system. We consider an incoming 
wave packet approaching the origin from x > and ask for the probability of crossing during 
a given time interval. We do this in two parts: first in this section, using a large time interval 
[0, r] and second in the next section, using a set of intervals of arbitrary size interval. 

A. Class Operators and Probabilities 

We consider first the following simple question. What is the probability of crossing or not 
crossing during the time interval [0, r]? The class operators for not crossing and crossing 
are 




(6.1) 




(6.2) 



and they satisfy 



e -iH r 



(6.3) 



We are interested in the probabilities for not crossing and crossing, 




(6.4) 
(6.5) 




and the off-diagonal term of the decoherence functional, 



D, 



Tr (C ncP Ct) 
Tr(C nc pe iHoT ) -p, 



nc 



(6.6) 



These quantities obey the relation 



Pnc + Pc + -Dc,nc + D* 



c,nc 



1 



(6.7) 
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We look for situations where there is decoherence, 

D c>nc = (6.8) 

(which is usually only approximate), in which case the probabilities then sum to 1, 

Pc(r)+ Pnc (r) = 1 (6.9) 

It is useful to relate some of these expression to the standard expressions for arrival time 
Il(t) defined in Eqs. fll.5p . (jl.6p (or Eg. (15. II) . To do this, note that p nc is in fact the same as 
the survival probability, N(t) defined in Eq. (ll.5[) . and that p nc obeys the trivial identity 

Pnc(r) = 1 + jTd* ^ (6.10) 

since p n c(0) = 1. It follows that 

p nc (r) = l- [ dtU(t) (6.11) 
Jo 

When there is decoherence, Eg. (16.91) holds and we may deduce that 

p c (r) = [ dt II(t) (6.12) 



Hence the decoherent histories analysis is compatible with the standard result, but only 
when there is decoherence. 



B. Calculation of the Decoherence Functional 



We now give two methods for checking for the decoherence of histories. The first involves 
expressing the probabilities and decoherence functional in terms of the transmitted and 
reflected waves defined in Eq. (I3.12[) . which implies that 

C n <#> = O(-xMr) + 0(X) (|We/> + (6-13) 
0#> = 6(-x) (|^/> - \lPtr)) ~ 9(xMref) (6.14) 

The probabilities and decoherence functional are therefore given by 

Pnc = (4>tr\iltr) + (Aef\Aef) + (Aef\^f) + (4>f\4>ref) + (ilf\0(x)\ijf) (6.15) 
Pc = {^tr\Ar) + (VVe/|VVe/) ~ {AM f) ~ (tpfl^tr) + (lpf\9{-x)\lj)f) (6.16) 
D C:nc = ('0t7#/) _ (ArWtr) ~ (AeflAef) ~ (^flAef) (6-17) 
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(Here for notational convenience we assume that the definition of \iptr) includes 9{—x) and 
that of |VVe/) includes 9(x), but the definition of the freely evolving part \ipf) does not 
include a 9 function). 

The magnitude of the off-diagonal term in the decoherence functional may be estimated 
from the explicit solution for the scattering states, Eqs. fl4.6p . (14. 8p . If there is substantial 
reflection, it is easily seen that the decoherence functional will not be small. So the inter- 
esting regime is the one explored in previous sections, namely Vo -C E (where E is a typical 
energy scale). In this regime the reflected wave functions are of order Vq/E. Furthermore, 
one can see from Eq.f l4.6p that the difference between \%p T ) and \iptr) are of order Vq/E. 
Therefore, the off-diagonal terms and the probability for not crossing are of order Vq/E, and 
the probability for crossing is of order 1, up to corrections of order Vq/E. Hence there is 
decoherence of histories in the regime Vq -C E. 

There is a second method of demonstrating decoherence which gives a different picture and 
will be useful later. Following the general pattern described in Eqs. ( 11.191) . ( 11.20p . consider 
the quantity, 

q nc (r) = Tr (C nc pe tH ° T ) (6.18) 
From Eq. ( 16.61) . we see that the decoherence functional may be written, 

D C:nc = g nc (r) - p nc (r) (6.19) 

This means that q nc = p nc when there is decoherence. Or to put it the other way round, 
decoherence of histories may be checked by comparing q nc with p nc and this is what we now 
do. Recall that p nc is given by Eqs. ( 16.101) . ( 16.111) (which hold in the absence of decoherence). 
We may write q nc in a similar form: 

gnc(r) = 1 + / dt-^ (6.20) 



The integrand is similar to H(t) defined in Eq. (l5.1l) . so we define 

dt 



U(t) = (6.21) 



We now have that 

qnc(r) = 1-1 dt fi(i) (6.22) 
Jo 

It then follows that the decoherence functional is 

D c , nc = J dt(u{t)-U{t)) (6.23) 
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To compute the decoherence functional we need to calculate U(t), which is given by 

n(t) = V (V>| exp(tH t) 9(-x) exp(-iH t-V 6(-x)t) \ip) (6.24) 

This is almost the same as n(t) except that the exponential on the left involves only H 
and not the complex potential (and also an overall factor of 2). We therefore follow the 
calculation of Il(t) in Section 5 with small modifications. With a little care, one may see 
that the final result is the same as that for IT(t), Eq. fl5.13l) . except that 2V is replaced with 
Vo, that is, 

f[(t) = V I ds e~ Vo{T - s) J{s) (6.25) 
Jo 

This result holds for timescales greater than 1/Vq and under the semi classical approximation 
Eq. ( 14.101) (which required E ^> Vq so is equivalent to the requirement of negligible reflection 
encountered above). Finally, a calculation similar to that of Eqs. ( 12.101 - ( 12. 131) implies that 

dt U(t) ta f dt J(t) (6.26) 
Jo 

as long as V r 3> 1. Since this result is independent of V , Il(t) will satisfy the same relation. 
We thus deduce that 

D c , nc w (6.27) 
hence there is decoherence, under the above conditions. 

VII. DECOHERENT HISTORIES ANALYSIS FOR AN ARBITRARY SET OF 
TIME INTERVALS 

We now turn to the more complicated question of much more refined histories, that may 
cross the origin at any one of a large number of times, during the time interval [0, r]. This 
corresponds more directly to the standard crossing probability, U(t)dt, the probability that 
the particle crosses during an infinitesimal time interval [t,t + dt]. 

A. Class Operators 

We have defined class operators Eqs. (16. II) . (16. 2\\ describing crossing or not crossing during 
a time interval [0, r]. We now split this time interval into n equal parts of size e, so r = ne 
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and we seek class operators describing crossing or not crossing during any one of the n 
intervals. We first note that 

e~ lHoe = C nc (e) + C c (e) (7.1) 

where C nc (e) and C c (e) are defined as in Eqs. fl6.ll) . f)6.2p except that here they are for a time 
interval [0, e]. We now use this to decompose e~ tH ° T into the desired class operators. We 
have 



— I ° / 

n-1 



= (e"^) n (C nc (e) + C c (e)) 
= (e-^)"" 1 C nc (e) + e^^^e) (7.2) 

Repeating the same steps on the first term, this yields, 
Repeating more times eventually yields, 

n— 1 

e~ lHoT = C nc (r) + Yl e- lHQ ^ k+1 ^C c (e)C nc (ke) (7.4) 

fc=0 

From this expression, we see that the class operator for crossing x = for the first time 
during the time interval [ke, (k + l)e] is given by the summand of the second term, 

C c ((k + l)e, ke) = e- lHo ^ k+1 ^C c (e)C nc (ke) (7.5) 

We will not in fact work with the class operator Eq. fl7.5p . since a more useful similar 
but alternative expression can also be found. Taking the continuum limit of Eq. fl7.4p (and 
inserting the explicit expression for C nc ), we obtain, 

e -iH r = e -iH r-Vr + f df e -iH (r-t)y e ~iH t-Vt ^ 

Jo 

(where, recall, V = Vq8(— x)). This indicates that the class operator for first crossing during 
the infinitesimal time interval [t, t + dt], is 

C c (t) = e -iHo(r-t)y e -iH t-Vt ( 77 ) 

We do not, however, expect histories characterized by such precise crossing time to be 
decoherent, so it is natural to consider coarser-grained class operators, 

C k c = [ L+1 dtC c (t) (7.8) 
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which represents crossing during one of the N time intervals [tk,tk+i] of size A, where 
tk = kA with k — 0, 1 • • • N — 1 and r = iVA. The complete set of class operators C a for 
crossing and not crossing is the set of iV + 1 operators 

C a = {C nc ,C k c } (7.9) 

and Eq. fl7.6p implies that they satisfy 

N-l 

e- m ° T = C nc +Y,C k c (7.10) 

k=0 

To check for decoherence of histories we need to calculate two types of decoherence functional 

D W = Tr (C c V(Cf )t) (7.11) 
D k>nc = Tr (C k c p{C nc y) (7.12) 

and this will be carried out below. 



B. An Important Simplification of the Class Operator 

There turns out to be a very useful simplification in the class operator Eq. (17.71) . Consider 
the amplitude 

(x\e iH ° T C c (t)\^) = V (x\e- iHot 6(-x)e~ iHot - vt \^} (7.13) 

for any x. The right-hand side is very similar to Eq. fl5.ip . except that there is no complex 
potential in one of the exponential terms and also the "final" state is \x) not \ip). (And 
there is also an overall factor of 2 different). Despite these differences, we may once again 
make use of the details of the calculation of Section 5, and we deduce from the analogous 
result Eq. fl5TT3T) . that 

(x\e iB ° T C e (t)W) = Vo f ds e~ Vo ^ ^(x\ (p 5{x s ) + 5(x s ) p) |V) (7.14) 

Jo 2m 

Like the derivation of Eq. (15.131) . this is valid under the conditions that all energy scales 
are much greater than Vq and all time scales much greater than 1/Vq. Now we integrate 
this over time to obtain the coarse-grained crossing time operator, Eq. fl7.8p . and again use 
approximations of the form Eqs. fl2.10l) - fl2.13p (again using the assumption of timescales much 
greater than 1/Vo), to yield, the remarkably simple and appealing form, 

^H r C k f h+1 dt (Zl) Qgfaj + 5 ^ t)f)) (7J5) 
Jt k 2m 
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This may also be written even more simply 



e iH r C k = Q ^ tk)) _ 0(£( ffc+1 )) 



(7.16) 



C. Probabilities for Crossing 

The above expressions for the crossing time class operator are the most important results 
of the paper and provide an immediate connection to the standard expression for the arrival 
time distribution. Supposing for the moment that there is decoherence of histories, we may 
assign probabilities to the histories. The probability for crossing during the time interval 



However, as noted in Eqs. fll.19l) . (11.201) when there is decoherence of histories, this expression 
for the probabilities for histories is equal to the simpler expression 



which is precisely the standard result! The expression for the probability q{t k ,t k+ i) is not 
positive in general (although is real in this it happens), but when there is deco- 

herence, it is equal to tfc+i) , which is positive. Hence the decoherent histories result 
coincides with the standard result under the somewhat special conditions of decoherence of 
histories. 

D. Decoherence of Histories and the Backflow Problem 

There is an interesting connection between decoherence of histories and backflow. To 
see this, consider the following simple case. We consider histories which either cross or do 
not cross the origin during the time interval [ti,^]- So the crossing and not crossing class 
operators are C and 1 — C, where 



Ffe j tk+i] is 




(7.17) 



q(t k ,t k+1 ) = Tr(C c V"° r ) 




(7.18) 




(7.19) 
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where we have adopted the notation Xk = x(tk) (and for convenience we have dropped the 
exponential factor which is just a matter of definition and drops out of all expression of 
interest). The decoherence functional is 

D = (C(l-C)) 

= (C) - (C 2 ) (7.20) 

This may also be written 

D = -((9(- Xl )9(x 2 ) + e(x 2 )9(-x x ))) (7.21) 

a form we will use below to check decoherence. When there is decoherence, D = and the 
probability for crossing is 

p(h,t 2 ) = (C 2 ) = (C) (7.22) 

As noted above, (C) is the standard result, Eq. fl7.18n . for the probability of crossing. 

There is an interesting connection here between backflow and decoherence. If there is 
decoherence, D is zero so (C) must cancel (C 2 ) in Eq.f l7.20p . which means that (C) > 0, so 
there is no backflow. Or we may make a logically equivalent statement: if there is backflow, 
(C) < 0, then there cannot be decoherence, since \D\ is then greater than the probability 
(C 2 ). Hence, states with backflow do not permit decoherence of histories. (Note that absence 
of backflow, (C) > 0, is not itself enough to guarantee decoherence - the stronger condition 
D = must be satisfied). 

This is an important result. The quantity (C) is regarded as the "standard" result for 
crossing time probability and its possible negativity is disturbing. Here, the decoherent 
histories approach sheds new light on this issue. In the decoherence histories apporach, the 
true probability for crossing is the manifestly positive quantity (C 2 ) and this is equal to 
(C) only when there is decoherence. In particular, when there is significant backflow, there 
cannot be decoherence, so probabilities cannot be assigned and (C 2 ) is not equal to (C). 

E. The Decoherence Conditions 

The crossing probabilities described above are only valid when all components of the 
decoherence functional, Eqs. fl7.lll) . (I7.12p . are zero. We therefore address the issue of finding 
those states for which there is negligible decoherence. 
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We consider first the simpler case, of the decoherence functional Eq. (17.12|) . The non- 
crossing class operator C nc is given in general by Eq. fl6.ip . However, it simplifies considerably 
in the approximations used to derive Eq. fl7.16p . which we adopt here. In particular, Eq. fl7.10j ) 
with Eq. (I7.16P imply that 

e lHoT = C nc + e iHoT [9{x) - 6{x{t))\ (7.23) 

Since we are interested only in initial states with support entirely in x > 0, we have 8(x) = 
which means that effectively, 

C nc (r) « 0{x)e- iHoT (7.24) 

The decoherence functional of interest is then 

D k , nc = (if>\6(x(T)) [6(x k ) - 9(x k+1 )} |V) (7.25) 

This is conveniently rewritten, 

D k , nc = W(x(t)) [6{-x k+l ) - e{-x k )\ |V) (7.26) 

We will take r to be very large and it is pretty clear that this object will be approximately 
zero, since we expect all the initial state to end up in x < at large times. However, we will 
see this below in more detail. 

The more important decoherence condition is that D kk i defined in Eq. fl7.lll) vanishes, so 
we now focus on that. We write the class operator Eq. fl7.16p for crossing during the time 
interval [t k , t k+ i] as 

C* = e- iHoT (6(x k )-6(x k+1 )) (7.27) 

For an initial state \ip), the quantity C^\ip) is a quantum state representing the property 
of crossing of the origin in the time interval [t k ,tk+i]- The decoherence condition D kk * = 
is simply the condition that the "crossing states" C^\ip) for different time intervals have 
negligible interference. The states C^\ip) consist of an initial state which has been localized 
to a range of time at x — 0. This is closely related to the interesting question of diffraction 



in time 



37j and this connection will be explored in more detail elsewhere 
The decoherence functional is given by 

D kj = ((e(x k )-e(x k+1 ))(e(x j )-e(x j+1 ))) 

= ((6(x k ) - 9(x k+1 )) - 0(- Xj ))) (7.28) 
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where without loss of generality we take tj + \ < tk- It is a sum of terms each of the form, 

d kj = (e{-x k )B{xj)) (7.29) 

where tk < tj. Note that 

14/ <d 2 m = (e(-x k )9(xj)e(-x k )) (7.30) 

The key thing is that d 2 m has the form of a probability - it is the probability to find the 
particle in x < at t k and then in x > at tj. Semiclassical expectations suggest that this 
is small in general for the states considered here, which are left-moving wave packets, and 
indeed the classical limit of this probability is zero. So this is a useful object to calculate 
in terms of checking decoherence. (Although it may not be small for states with backflow). 
Note that it also implies that the other parts of the decoherence functional, Eq. (17.261) will 
also be small. In detailed calculations, the upper bound Eq. ( 17.30 1 ) must be compared with 
the probabilities, as in Eq. (1 1.1 81) . 



F. Checking the Decoherence Condition for Wave Packets 

We now consider the particular case of an initial state consisting of a wave packet 

= ( 27 raV /4 6XP (- ^J^ + mX ) (7 - 31) 
where go > and po < 0. We first consider a heuristic analysis of decoherence. In the 
simplest case, the wave packet crosses the origin almost entirely during the time interval 
[tk, tk+i] (of size A) for some k, without any substantial overlap with any other time intervals. 
(See Figure 1.) This means that 

Cf « for k'^k (7.32) 

and it follows that Df.y ~ 0. The key time scales here are the classical arrival time for the 
centre of the packet, 



and the Zeno time 



bo I 



ma 

u = w (7 ' 34) 
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(which is also approximately equal to 1/(AHq)). The Zeno time is the time taken for the 
wave packet to move a distance equal to it spatial width a, or equivalently, it is the size of 
the packet's "temporal imprint" at the origin. Therefore, the above approximations work if, 
firstly, 

t g < A (7.35) 

and secondly, if the classical arrival time t a lies inside the interval [tfe, tfc+i] and is at least 
one or two Zeno times away from the boundaries. 

It is easy to see that similar conclusions hold for superpositions of initial states of the 
form Eq. (j7.3ip as long as they are approximately orthogonal. Loosely, this is because under 
the above conditions, the class operators do not disturb the states and the only non-zero 
components of the off-diagonal terms of the decoherence functional will be proportional to 
the overlap of pairs of initial wave packets, so will be approximately zero. (See Figure 3.) 
More general, non-orthogonal superpositions may, however, produce backflow, so there may 
be no decoherence. 



G. Checking Decoherence for a Detailed Model 

Decoherence starts to become lost as the the size A of the time intervals [tft, tfe+i] is 
reduced to close to the Zeno time. This is because the wave packet will split into parts that 
cross during different time intervals and the effect of diffraction in time mentioned above 37] 
will cause these different parts to be non-orthogonal. These effects will be explored in more 
detail elsewhere 38|. Here, we give a more detailed calculation to check for decoherence. 

For simplicity we work with the simple case considered in Subsection 7D above. We take 
the initial state to be the wave packet Eq. (l7.3ip and we note that the decoherence functional 
Eg. (17211) satisfies 

\D\ 2 < 2d 2 m (7.36) 

with d 2 m given by Eq. (17.301) with k — 1, j — 2. We need some probabilities to compare this 
with. We have that 

\D\ 2 < (C 2 )((l-C) 2 ) (7.37) 

The interesting case is that in which the crossing probability (C 2 ) is somewhat less than 1, 
less than about 1/2, say, in which case the non-crossing probability ((1 — C) 2 ) will be of 
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order 1. It is therefore sufficient to compare d 2 m with (C 2 ). Now note that 

(c 2 ) = ((0(xO - e{x 2 )f) 

= (0(zi)0(-z 2 )0(^i)) + (e(zi)0(-x 2 )0(-£i)) + (0(z 2 )0(-£i)) (7.38) 

This means that (C 2 ) is in fact equal to the probability 

P12 = (fl(xiM-x 2 )fl(xi)> (7.39) 

up to terms which vanish when D — 0. This is useful since it is now identical in form to the 
expression for d 2 m and our goal is to show that 

d 2 m < Pu (7.40) 



It is useful to work in the Wigner representation [39], defined, for a state p(x,y) by 

W(p, q) = ±- I ' di e'^ p(q + U,q- U) (7.41) 
In J 11 

The probabilities pi 2 and d 2 m are then given by 



pi 2 = 2tt / rfpdg Wi a (p,g) W (p,g,ti) (7.42) 
d 2 m = 2n J dpdq W D {p,q) W Q {p,q,h) (7.43) 
Here, Wo(p, g, ti) is the Wigner function of the initial state, evolved in time to t\, 

2a 2 

The objects Wp and W^o are the Wigner transforms of the ^-function combinations appearing 
in Eqs.f l7.30j) . (17.391) and are given by 



W (p, q, h) = — exp I — 2a (p - p ) ] (7.44) 



W P (p,q) = ^9{q) f dy S -^- (7.45) 

2 ^ Ju(p,q) y 

W D (p,q) = ±6(-q) r dy S -^ (7.46) 
2vH J uM y 



u(p,q) = 2q[p+ . mq . I (7.47) 



where 

(*2-tl). 

We see that the only difference between the expressions for pi 2 and d 2 m is in the sign in the 
^-functions. 
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The integral 



f(u) = I dy S -^ (7.48) 

y 



may be expressed in terms of the Sine integral function Si(x), 

f(u) = \- Si(u) (7.49) 

but its properties are not hard to see directly. For large negative u, f(u) ~ ir, at u = 0, 
/(0) = 7r/2, and for large positive u, f(u) goes to zero, oscillating around 1/u. (See Figure 
4.) 

We now compare the size of pi 2 and d^.We assume that the wave packet is spatially broad, 
so a is large and the Wigner function Eq. fl7.44p is therefore concentrated strongly about 
P — Po < 0. We therefore integrate out p and set p = p throughout. The most important 
case to check is that in which the wave packet is reasonably evenly divided between x > 
and x < at time t\ so that both p 12 and d 2 m have a chance of being reasonably large. This 
means that go + Po^i/ m should be close to zero (to within a few widths a), so for simplicity 
we take it to be exactly zero. 

With these simplifications, we have 

Here, since q < and po < 0, we have u(po,q) > 0. We can evaluate this expression by 
examining the comparative effects of f(u) and the exponential term. From the plot of f(u) 
(see Figure 4), we see that it drops to zero at u = w (which is of order 1) and oscillates 
rapidly around zero for m > «o, so we expect the integral to be dominated by values of q for 
which < u < u . The value u = u corresponds to q = q where 



= bo I A 
2m 



1+ Uo 



1/2 

- 1 I (7.51) 



E A_ 

where E = p\j1m. In the complex potential calculations, we have assumed that E ^> Vq 
and we also assumed that all timescales are much greater than 1/Vo, and these together 
imply that E A ^> 1. We may therefore expand the square root to leading order and obtain 

* K ~2R (7 ' 52) 

We have assumed that the wave packet is sufficiently broad that apo 3> 1, and this means 
that 

M < (7-53) 
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This in turn means that f(u) is significantly different from zero only in the range \q\ -C a, and 
most importantly, in this range, the exponential term in Eq. (I7.50j) is approximately constant. 
We may therefore evaluate Eq. (I7.5UI) by ignoring the exponential term, integrating from 
to go an d approximating f(u) as 

f(u)^^-u + 0(u 3 ) (7.54) 

We thus obtain the simple result, 

2 1 1 

In the expression for p 12 , there is a key difference in that q > which means that u(po, q) 
can be positive or negative. Introducing 

IPol a a 

q z = A = a— (7.56) 

m t z 

(where t z is the Zeno time) we see that u < for q < q z and u > for q > q z . We therefore 
have 



P* = (2^2)1/2 [ d * ex p("2^) 



+ / Am/9 / d ? ex P f-7TT ) f( u (Po,q)) (<•••""•") 



(2ttV 2 ) 1 /2 ^ ^ ^ 2 (T 2 

Here, f(u(po,q)) ~ vr in the first term, differing from this value only in a region of size 1/po 
close to q — 0. In the second term f(u) will tend to be small except for a small region of 
size 1/po around the origin. 

If q z ^> cr, (that is, A ^> t z ), the second term in p 12 is exponentially suppressed and in 
the first term the integration range is effectively to oo, so we obtain 

Pn « ^ (7.58) 

This is the expected result, since under the above assumptions on the wave packet, half of 
it will cross x — if the time interval is sufficiently large. Clearly p 12 3> in this case so 
there is decoherence. 

If q z a, the first term in p\ 2 is of order q z /a and the second of order l/( |po ? the 
same order of magnitude as d? m . Hence in this case we have decoherence if 

A 1 , 

T » rr 7.59 
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Since the right-hand side is already <C 1, this is easily satisfied even for time intervals whose 
size is A is smaller than the Zeno time. In fact this condition is equivalent to the condition 
E A 3> 1, which is satisfied by the assumptions of the complex potential model, as stated 
above. 

In brief, we therefore get decoherence of histories for a single wave packet under a wide 
variety of circumstances. 

VIII. SUMMARY AND CONCLUSIONS 

This paper was initially motivated by a desire to analyse the arrival time problem using 
the decoherent histories approach to quantum theory. But along the way we have reconsid- 
ered and derived a number of other useful related results. 

We considered the arrival time problem using a complex potential to kill paths entering 
x < 0. In Section 2 we gave a classical analysis of the problem. We derived a result of the 
expected form exposing the resolution function as an essentially classical effect summarizing 
the role of the complex potential. We also showed that coarse graining over time scales 
much greater than 1/Vb produces a formula for the arrival time of expect form and which is 
independent of the complex potential. This is an important result for the rest of the paper. 

After reviewing the path decomposition expansion in Section 3, we used it in Section 4 to 
derive the standard results for stationary scattering states. This is simpler and more direct 
than the usual method, involving matching eigenfunctions. (In a sense this cumbersome 
machinery is already encoded in the PDX). In Section 5, we used the PDX to rederive the 
standard form of the arrival time distribution with a complex potential, in the limit of weak 
potential. The form of this calculation turned out to be useful for the subsequent work on 
the decoherent histories approach. 

In Section 6, we considered the decoherent histories analysis for the simple case of a 
particle crossing or not crossing x = during a large time interval [0, r]. We found the 
simple and expected result that the histories are decoherent as long as reflection by the 
complex potential is negligible. The resultant probabilities are consistent with the standard 
result for the arrival time. 

The main part of the decoherent histories analysis was given in Section 7, where we 
first derived the class operators described crossing x = for an arbitrary set of small time 
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intervals. Here we obtained our most important result: the crossing class operator Eq. fl7.15p 
for timescales much greater than 1/Vq. This form of the class operator gives an immediate 
connection with the standard result for probabilities when there is decoherence. Indeed, one 
may have guessed the form of the class operator from the standard form of the probabilities, 
and this is pursued in another paper 30| . However, it is also gratifying that it can be derived 
in some detail using the complex potential approach used here. 

To assign probabilities, the decoherence functional must be diagonal and we considered 
this condition. We found a variety of states for which there is decoherence, under certain 
more detailed conditions, which we discussed. 

We also noted an interesting and important relationship between decoherence and back- 
flow: If there is decoherence, the probabilities for crossing must be positive so there cannot 
be any backflow. If there is no decoherence, the integrated current may still be positive, but 
one can say that if there is backflow there will definitely be no decoherence. This means 
that the decoherent histories approach brings something genuinely new to the arrival time 
problem: it establishes the conditions under which probabilities can be assigned and in 
particular forbids the assignment of probabilities in cases where there is backflow. 
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APPENDIX A: SOME PROPERTIES OF THE CURRENT 

We have derived the expression 

p(0,T) = [ dt J(t) (Al) 
Jo 

as the approximate probability for crossing the origin during the time interval [0, T], where 
J(t) is the usual quantum-mechanical current. The current itself is not necessarily positive 
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due to backflow. Here we explore the possibility that averaging it over time might improve 
the situation. On the one hand, the results of Bracken and Melloy [3] show that there is 
always some state for which p(0, T) defined above is negative, for any T. On the other 
hand, for a given state, one might hope that p(0,T) will be positive for sufficiently large T. 
Here we give a brief argument for this, which also makes contact with the negativity of the 
Wigner function. 

The current can be written in terms of the Wigner function W(p,x) as 

J(t) = - [ dp — W(p,0,t). (A2) 
J m 

The Wigner function evolves freely according to W (p, x, t) = W(p, x—pt/m,0). We assume 
it has support only on negative momentum states, with average momentum p < and 
momentum width o~ p . 

Consider a time interval < t < T over which backflow occurs. It is clear that in order 
for this to occur the Wigner function must be negative for at least some of this interval. We 
can write, 



p(0,T) = - f dt f dp — W(p,-pt/m,0) 
Jo J m 

/r\p\T/m 
dp J dx W(p,x,0) (A3) 



So p(0, T) is given by the average of the Wigner function over a region of phase space. We 
now recall a standard property of the Wigner function which is that, broadly speaking, it 
will tend to be positive when averaged over a region of phase space of size greater than 
order 1 (in the units used here where H = 1). This region is of size of order \p \T/m in the 
x-direction but infinite in the p-direction. However, the Wigner function has momentum 
spread a p , so the effective size averaged over is a p p T/m which is approximately the same 
as AH T. This means that we expect that p(0, T) will be positive as long as 

T > ~ee (A4) 

Hence, as expected, the integrated current will be positive for T sufficiently large and the 
key timescale is the Zeno time. This heuristic argument will be revisited in more detail 
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elsewhere. 
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Figure 1. The quantum arrival time problem. We prepare an initial state localized entirely 
in x > and consisting entirely of negative momenta. What is the probability that the 
particle crosses the origin during the time interval [£1,^2]? 
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Figure 2. The path decomposition expansion. Any path from xq > at t = to a 
final point x\ < at t = T has a first crossing of x — at t\ and a last crossing at t^. 
The propagator from (x ,0) to (xi,T) may be decomposed into three parts: (A) restricted 
propagation entirely in x > 0, (B) free propagation starting and ending on x = 0, and (C) 
restricted propagation entirely in x < 0. The corresponding path decomposition expansion 
formula is given in Eq. (13.61) . 
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Figure 3. An initial state consisting of a superposition of wave packets may have significant 
crossings in at least two different time intervals. If the initial packets are orthogonal, and 
the time intervals are sufficiently large (greater than the Zeno time), the packets will remain 
orthogonal after passing through the time intervals and the corresponding histories will be 
decoherent. 
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Figure 4. The function f(u). It oscillates around zero for u 3> and oscillates around n 
for n < 0. As a function of q, f(u(po, q)) differs from or n only in a region of size l/\po\ 
around q = 0. 
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